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ABSTRACT 


A  new  type  of  digital  filtering  it  explained  without  ualng  the  convolution  theorem  used  in  electrical 
engineering.  The  icheme  makes  use  of  the  Numerical  Transform  Theorem  which  transforms  a 
pair  of  finite  spectra  (1.  e. ,  amplitude  and  phase  spectra)  in  the  frequency  domain  into  a  finite 
time  series.  This  time  series,  which  may  be  filtered  during  transformation,  has  precisely  the 
same  record  length  as  the  "original"  time  series  without  introducing  meaningless  transients.  The 
spectrum  filtering  technique,  in  contrast  to  the  convolution  filtering  technique,  has  an  additional 
advantage  that  it  does  not  require  a  separate  digital  filter,  Some  examples  are  given* 


SPECTRUM  FILTERING  BY  MEANS  OF 
NUMERICAL  TRANSFORM  THEOREM 


».  INTRODUCTION 


For  many  year*,  Fourier  aerie*  and  traneform*  have  been  u*ed  in  (cienee  and  engineering.  The 
primary  reaeon  lie*  in  it*  simplicity  of  representation*  which  allow  ua  to  analyze  a  function  in  a 
different  form. 

Foi  purposes  of  numerical  analyses,  series  truncation  becomes  inevitable,  and  expressions  for  a 
■  et  of  coefficients  are  found.  There  1*  an  important  theorem,  however,  which  has  escaped  the 
attention  of  previous  investigator*.  This  theorem,  which  may  be  called  "Numerical  Transform 
Theorem",  transforms  finite  time  series  to  finite  frequency  harmonics. 

The  history  of  frequency  filtering  is  probably  as  old  as  the  history  of  electronic  engineering.  So 
far,  digital  filtering  it  accomplished  through  convolution  theorem  which  states  that  a  multiplies* 
tlon  in  a  frequency  domain  corresponds  to  a  convolution  in  the  time  domain.  Theoretically,  the 
convolution  filter  has  infinite  length,  and  a  truncation  of  this  filter  at  a  certain  tolerable  level  i* 
necessitated  in  actual  applications.  The  resulting  filtered  record  contains  useless  tails  caused  by 
transients  in  the  convolution  process. 

This  article  intends  to  show  a  different  way  of  spectrum  filtering  by  which  a  time  series  can  be 
filtered  without  generating  convolution  transients.  A  terminology  "Spectrum  Filtering"  it  adopted 
to  differentiate  this  technique  from  the  conventional  convolution  filtering  technique. 

This  work  was  performed  under  Project  VT/4051,  Contract  AF  33(657)-12145.  The  project  is 
under  the  technical  direction  of  the  Air  Force  Technical  Application*  Center  (AFTAC)  and  the 
overall  direction  of  the  Advanced  Research  Projects  Agency  (ARPA). 


2.  NUMERICAL,  TRANSFORM  THEOREM 


Since  the  time  Fourier  made  use  of  the  trigonometric  series  beating  his  name  for  the  Integration 
of  the  equations  of  thermal  conduction,  the  technique  has  become  one  of  the  most  basic  tools  in 
scientific  and  engineering  analyses.  The  method  involve*  an  expression  of  a  continuous  function 
of  length  2T  by  an  infinite  series, 


where  coefficients  a'n  and  b'n  can  be  determined  by 


-i-  |ZT  f(t)  cos  dt 

I  ft 


l  nrt 


(l) 


(2) 


Although  these  expressions  are  exact  by  virtue  of  the  orthogonality  relationship  of  trigonometric 
functions,  it  takes  an  infinite  number  of  *'n  and  b'n  to  have  an  exact  expression  of  the  original  f(t). 


There  1*  an  important  theorem  in  the  theory  of  tranaformii  which  la  uiually  referred  to  ae 
Fourier  tranaform  theorem  (Lee*  1961,  p.  33).  The  mathematical  exer*aa<nr>  of  this  theorem  is 
given  by 

f(t)  •  j  df  J  Hi)  e  ‘J2rT£  (3) 

where  the  frequency  la  given  in  cydea  per  aeeond. 

This  expresalon  holds  for  continuous  function*  extending  from  *“  to  *  in  both  time  and  frequency 
domains.  The  approximation  of  the  above  expression  in  numerical  form  takes  infinite  terms,  and 
it  is  difficult  to  recover  the  original  function  by  using  the  equation  (3). 

In  my  previous  preaantationl  at  St.  Louis,  Missouri,  1  indicated  a  theorem  called  Numerical 
Transform  Theorem  (Huang,  196$).  Borrowing  the  notation  of  Z-transform  (Robinson  and  Treltel, 
1964),  this  theorem  can  be  expressed  by 


where  -j-fcp  (3) 

Zee 

In  the  above  expressions,  J  is  an  imaginary  unit  and  m  is  the  sample  number.  The  expression  (4) 
is  also  an  exact  relationship  which  transforms  a  finite  number  of  digitized  data  into  a  finite 
domain.  It  indicates  that  the  exact  recovery  of  the  original  digitized  data  ie  possible  without 
consuming  an  infinite  number  of  coefficients. 


3.  SPECTRUM  FILTERING 

By  decomposing  equation  (3)  into  a  transform  pair,  we  have 


F(f)  - 

r  , 

-JStff 

/  f(f>  « 

df 

f(t)  . 

(*  F(f)  • 

^£t  df 

and 


Denoting  a  filter  function  by  H(f),  the  impulse  response  of  this  function  is  given  by 


(6) 

V) 


h(t) 


•£ 


)2nft 

H(f)  a  df 


(8) 


^Presented  at  the  Seismologicel  Society  of  America  meeting,  St.  Louie,  Missouri,  on 
14  April  1963. 
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Figure  la.  Analog  d 
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Figure  lb.  Digitized  deep-^eil 


For  a  filter  length  of  pAt,  the  fundamental  frequency  is  given  by 


Any  frequency  f  in  terms  of  Af  can  be  expressed  by 


The  mathematical  representation  of  a  low-pass  filter  is 


(15) 


(16’ 


(1,  for  n  =  0,  1, 2, .  .  .  N(N  *  q) 

(17) 

0,  for  n  >  N 

The  impulse  response  of  this  filter  is 
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ta 


-ni 


(18) 


Figure  2  shows  the  analog  version  of  the  low-pass  filtered  record  (24  dB/octave  at  1.5  cps)  of  the 
original  figure  la.  An  impulse  response  of  the  low-pass  filter  (cut-off  at  1.5  cps)  is  shown  in 
figure  3a  and  the  digitally  filtered  result  using  the  conventional  convolution  technique  is  shown  in 
figure  3b. 


Figure  4  is  the  corresponding  low-pass  filtered  record  by  a  direct  usage  of  the  numerical 
transform  theorem. 

b.  High-pass  filter  — 

For  a  high-pass  filter,  the  filtering  function  becomes 


0,  for  n  “  0, 1,  2, .  .  .  N  (N  <  q) 

Hn  = 

1,  for  n  >  N 

The  impulse  response  of  thiB  filter  is 


(19) 


ht  a  Z'ni  (20) 

Figure  5  is  the  analog  high-pass  filtered  result.  The  ...ipulse  response  of  the  high-pass  filter 
together  with  the  convolved  result  is  shown  in  figures  6a  and  6b.  Figure  7  is  the  filtered  output 
by  using  the  numerical  transform  theorem. 


c.  Band-pass  filter 

The  mathematical  expression  for  a  band-paBS  filter  is  given  by 


H 

n 


0,  for  n  e  0,  1,  2,  .  .  .  Nj 

1,  for  n  =  +  1,  N  j  +  2,  .  .  .  W2 

0 ,  fo  r  n  >  N  2 


(21) 
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Figure  7.  High-pass  filtered  Apache  record  using  the  numerical 


ui  this  UJtcr  is  accordingly 


h.  « 

l 


z-ni 


Figures  8  through  10  show  different  results  of  the  filtered  outputs. 


(22) 


As  a  comparison  between  the  analog  results  and  the  reconstructed  results,  figure  11  and  figure  12 
are  displaye  It  is  of  interest  to  note  that  except  for  the  dc  bias,  figure  12  is  coincident  with  the 
figure  lb  on  which  the  spectrum  filterings  were  performed.  Other  types  of  filteringi  e.g. ,  notch 
filter,  comb  filter,  etc.,  can  easily  be  arranged  using  the  numerical  transform  theorem. 


5.  CONCLUSIONS 


By  virtue  of  the  numerical  transform  theorem,  it  is  possible  to  filter  digitized  records  in  any 
desired  fashion  in  the  frequency  domain.  No  such  useless  transients  are  incurred  in  this  process 
as  in  the  case  of  filtering  through  convolutions. 

When  the  filtered  function  using  the  numerical  transform  theorem  is  spectral  analyzed  using  the 
Fourier  transform,  there  will  be  small  side  lobes  created  between  the  multiples  of  the  original 
fundamental  frequency.  This  situation  can  be  remedied  by  taking  a  sufficiently  long  record 
length  Initially,  The  increased  frequency  resolution  will  thus  contribute  to  a  reduction  in  the 
undesirable  side  lobes.  This  is  the  main  trade-off  between  the  analog  and  the  digital  filterings. 

A  comparison  between  figure  11  and  figure  12  indicates  that  digital  filtering  is  closer  to  a  perfect 
filter  than  analog,  and  the  application  of  the  numerical  transform  theorem  will  yeild  filtered 
records  without  transients.  The  transient  effects  are  particularly  significant  for  records  of 
short  durations. 

Although  one -dimensional  transformation  is  shown  in  this  article,  the  extention  of  this  filtering 
technique  can  be  accomplished  to  three-dimensional  array  problems  where  both  time  and  space 
parameters  must  be  taken  Into  consideration.  It  must  be  remembered  that  the  above  technique 
requires  spectral  analysis  before  filtering,  although  a  separate  digital  filter  la  not  required  in 
the  procees. 
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